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ABSTRACT 

We combine two approaches to isolate the AGN luminosity at near-infrared wavelengths and relate the near-IR pure 
AGN luminosity to other tracers of the AGN. Using integral-field spectroscopic data of an archival sample of 51 local 
AGNs, we estimate the fraction of non-stellar light by comparing the nuclear equivalent width of the stellar 2.3 jum GO 
absorption feature with the intrinsic value for each galaxy. We compare this fraction to that derived from a spectral 
decomposition of the integrated light in the central arc second and find them to be consistent with each other. Using 
our estimates of the near-IR AGN light, we find a strong correlation with presumably isotropic AGN tracers. We show 
that a significant offset exists between type 1 and type 2 sources in the sense that type 1 sources are 7 (10) times 
brighter in the near-IR at logLACN = 42.5 (logLACN = 42.5). These offsets only becomes clear when treating infrared 
type 1 sources as type 1 AGNs. 

All AGNs have very red near-to-mid-IR dust colors. This, as well as the range of observed near-IR temperatures, can be 
explained with a simple model with only two free parameters: the obscuration to the hot dust and the ratio between the 
warm and hot dust areas. We find obscurations of = 5... 15 mag for infrared type 1 sources and = 15... 35 
mag for type 2 sources. The ratio of hot dust to warm dust areas of about 1000 is nicely consistent with the ratio of 
radii of the respective regions as found by infrared interferometry. 
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1. Introduction 

Most, if not all, Active Galactic Nuclei (AGNs) have nu¬ 
clear dust structures on parsec scales that are made re- 
spous ible for the line of sight obscuration in edge-on sys¬ 


tems (Antonucci 
radiation from t’ 
these “tori' 


1993) and the collimation of the ionizing 

IMl. While 


re AGN (e.g. Fischer et al 
were initially invoked to explain the differ¬ 


ences between edge-on and face-on systems, it is increas¬ 
ingly realized that they are not simple “donuts”. Instead, 
they are part of a complex and dynamical accre tion flow 
towards the central super-massive black hole (e.g. Norman 


fc Scoville|1988[ Wada fc Norman|2QQ2{|Vollmer et al.|2QQ8t 


Wada et al.|2Q09[|Schartmann et al 


Quataert||26lOt |Hopkins et al.||2012|). Multi-scale hydrody- 

namical simulations are successful in predicting observed 
properties such as t he scale-size and inferred mas s of the 
parsec-scale “torus” ( Schartmann et al.|2Q09 ,|2QlQ[) and ex¬ 


wavelengths (e.g. Prieto et al. 2010). While most tori are 


unresolved even in the highest resolution single-dish obser¬ 


vations (e.g. [Asmus et al.[2014 ), more than two dozen have 


now been spatially resolved using near- and mid-infrared 
interf erometry (e.g. [Kishimoto et al.||2011 Burtscher et al. 
2013) and their mia-infrared sizes have been measured to 


range between about 0.1 and 10 pc for Seyfert galaxies and 
quasars up to a bolometric luminosity of Lboi = 10^^ erg/s. 

In infrared observations, the most prominent difference 
between face-on (type 1) and edge-on (type 2) systems 
is seen in the Silicate feature at 9.7 /im which is often 
in emission i n type 1 and most ly in absorption in type 2 


obje cts (e.g. Honig et al. |2010 ). Models with continuous 
(e.g. IScharto ann et al.||2005|) and clumpy (e.g. jNenkova 
" ' ib 


'opkins 


et al.||2002[) dust distributions have been invoked to ex¬ 
plain tEe'^served differences between the various classes 
of AGNs and to relate them to torus properties such as 
inclination angle, number of clumps or opening angle (see 
HQnig|20T3 for a recent review). It is unclear, however, how 


plain the appearance of dust y disk structures as a dynami¬ 
cal consequence of accretion ( [Hopkins et al.|[2Q12 ). 

This region (we will keep calling it torus for the sake of 
simplicity) is very prominently seen by its thermal emission. 
It re-radiates a substantial fraction of the bolometric lumi¬ 
nosity of the AG N central engine in near- to mid-infrared 


accurately torus parameters can be retrieved from SED fits 
as the assumptions about the torus structure (e.g. power- 
law distribution in clump density) may not be correct. In¬ 
terferometric observations have found that essentially all 
infrared bright tori consist of two sp atially distinct compo- 
nents which are very different in size ( Burtscher et al.|2Q13 ) 
which seems to contradict the simple cloud distributions as- 
sumed i n models. Moreover , in a comparison between torus 
models, Feltre et al. (2012) found that most of the differ¬ 
ences in the SEDs between models arise not from the dust 
morphology but from the model assumptions. 
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On the other hand, all models predict significantly dif¬ 
ferent intensities from dust emission at infrared wavelengths 
between edge-on and face-on inclinations. In many stud¬ 


ies, however, no such differences were observ ed (Lutz et al. 
2004 Gandhi et al.||2009 Asmns et al.||2011 ), thereby chal- 


mc. 


ination-dependent AGN unifica- 


lenging the classic^ 
tion scheme. 

In this study we compile a large archival sample of high- 
resolution near-IR integral field spectroscopic data of local 
AGNs and determine the fraction of non-stellar light in the 
near-IR with two complementary methods. We describe the 
sample and the data analysis in Section and the decom¬ 
position of the nuclear light into AGN and stellar light in 
Section We then compare the near-IR AGN luminos¬ 
ity with other tracers of the AGN luminosity and discuss 
the luminosity relations in the context of obscuration and 
isotropy (Section |^. We summarize our findings in Section 


El 


2. Sample selection and data reduction 

We sel ect observations of all objects from the ES Q SIN- 
FONI (jEisenhauer et al.||2003| |Bonnet et al.||2004) archive 
that have an AGN classincation in SIMBAD. Additionally 
we impose a distance limit (< 60 Mpc) in order to well re¬ 
solve the nuclear region of the AGN. At the maximum dis¬ 
tance, our resolution is about 30 pc (300 pc) in the medium 
resolution cubes with (without) AO. We searched the ESO 
archive for SINEONI observations in the H+K or K grat¬ 
ings with any pixel scale and required only 10 minutes or 
more of time on target with seeing better than about 1 arc- 
secQln addition we add one source observed with Gemini 
NIES (NGC 4051), the reduced cube of which was kindly 
provided to us by R. A. Riffel. This resulted in a sample of 
51 active galaxies. The list of targets and observations can 
be found in Tab. [H 

Out of the 51 objects, 40 have been observed with AO (6 
with laser guide star) and 11 without. We detect non-stellar 
continuum associated with the AGN in the near-IR in 31 
sources (see Section^. Only five of these sources have been 
observed without AO. The median seeing of these cubes is 
0.9 arcsec. 

Eor the distance we adopt the median value of redshift- 
independent distance measurements from NED, where 
available, else the luminosity distance computed by NED 
using a ACDM cosmology with Hq = 73 km/s/Mpc, 
^matter = 0.27, f^vacuum = 0.73. Eor thosc galaxics where 
our X-ray reference (see below) quotes luminosities we 
adopt the distance used in the particular reference. 


^ To find these observations, we downloaded the headers of all 
SINFONI observations from the ESO archive and parsed the 
relevant data into an SQLite database. From that we constructed 
a list of unique coordinates which we parsed through SIMBAD 
to retrieve the object identifications and classifications based on 
these coordinates. From this list we filtered out local AGNs with 
sufficient data taken with appropriate settings, sufficient signal 
to noise and under good atmospheric conditions. 
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Data reduction was performed with spred, a custom 
packa ge developed at M PE for the analysis of SINFONI 
data ( Abuter et al.|2Q06 ). For th e sky subtrac tion, we used 
the routines mxcor and skysub ( Davies|[2QQ7 ) which corre¬ 
late reconstructed object and sky cubes spectrally and shift 
them such that their OH line wavelengths match before sub¬ 
tracting the scaled sky frame. This reduces sky subtraction 
residuals that otherwise occur because of a change in the 
wavelength scale between target and sky observations. This 
was followed by bad pixel and cosmic ray treatment using 
lac3d, a 3D L aplacian Edge Dete ction algorithm based on 
L.A.COSMIC ( van DokkumpQOl ). 

The typical calibration uncertainty is better than 0.1 
mag (or about 10 % in flux) as determined from the stan¬ 
dard deviation of aperture photometries of the individ¬ 
ual observations per source (before combining). The major 
source of uncertainty is variations in the conversion fac¬ 
tor within and between nights. Only for NGC 7743 is the 
calibration uncertain to about 0.5 mag because the obser¬ 
vations were performed under variable conditions. However, 
this source does not show any signs for AGN emission in 
the near-IR and is therefore not used in our quantitative 
analysis. 

2.1. Mid-IR, X-ray and [0 IV] samples 

In order to compare our near-IR derived AGN luminosities 
with other AGN indicators we collect X-ray and mid-IR 
data for all sources. Ground-based, high spatial resolution 
observations in the mid-infrared N band (8-13 pm) ex¬ 
ist for all objects with detected near-IR AGN light except 
for NGC 7496 and also for about half of the near-IR non¬ 
detections. The mid-IR data are mostly from VLT/VISIR 
observations whose high spatial resolution ensures mini¬ 
mum contamination by the host galaxy. Most of the near-IR 
detected ga laxies (24) also appear in the 70-month BAT all¬ 
sky survey (Baumgartner et al. 2013) in the very hard X-ray 
band 14-195 keV. This waveband is particularly helpful for 
AGN studies since it is essentially unaffected by obscura¬ 
tion except for the most Compton-thick objects. For six 
further objects, absorption-corrected 2-10 keV fluxes have 
been published which we use with the conversion factor 
Ai4-i 9 5keV ^ 5 X I/2-iokeV taken from Fig. 6 of IWinter 


et al. (2009). NGC 7496 is again the odd-one out which 
does not have a published hard X-ray flux although we de¬ 
tect non-stellar light in the near-IR. Conversely, only two 
near-IR non-detected sources are BAT-detected, but about 
half of the near-IR non-detections have published mid-IR 
fluxes. 

Additionally, we collected [OIV]A25.89/im observations 
for those objects with near-IR detections from the litera¬ 
ture. This mid-infrared line is only little affected by dust 
and due to its high ionization potential of 54.9 eV it is also 
essentially unaffected by starlight. It has been shown to be 
a goo d isotropic tracer of the AG N luminosity ([Rigby et al. 


lution of the 2.3 pm CO absorption feature with the 

nuclear mid-IR luminosity with the hard X-ray lu¬ 
minosity with and the [OIV]A25.89/im luminosity by 

^AGN • collection of multi-wavelength data is presented 

in Tab. O 

Errors for these luminosities are taken from the respec¬ 
tive references. The median errors on 
are 0.06, 0.02, 0.03 dex respectively; error bars are shown 
in the plots, but omitted in the Table for compactness. 


3. Stellar and AGN light decomposition, 
photometry and auxiliary data 

We decompose the stellar and AGN light in the near-IR by 

two complementary methods: 

1. Using the radial profile of the stellar GO absorption fea- 
ture and integrating in a nuclear aperture (Section |3.1| ) 

2. By fitting a stellar template and a blackbody to the 
spectrum extracted in the same nuclear aperture (Sec- 

Iq oh 


tion 3.2) 


3.1. Radial decomposition: The CO equivalent width radial 
profile 

First we smooth the cubes by a simple 5x5 pixel boxcar 
filter in order to improve the spectral signal to noise and 
increase the significance of the fit. This comes at the cost of 
spatial resolution, but since we are averaging over an aper¬ 
ture larger than this at the end, this does not impact our 
results. We then use the custom-built routine starfit^lto 
fit a stellar template to the CO(2,0) stellar absorption fea¬ 
ture at 2.29 pm at every spatial pixel of the smoothed cubes 
to obtain the stellar kinematics. We used a SINFONI spec¬ 
trum of the Ml star HD 176617. The specific star does not 
matter, however, since we not only normalize the continuum 
level, but subtract it, so that the continuum is effectively 
set to zero and the EW of the CO featur e has much less im- 
pact on the fit than in other methods; see|Engel et al. (2010) 
for a detailed discussion on this. With the known velocity 
field, we determine the equivalent width (E W) of this CO 
feature in the standard wavelength range (e.g. [Origlia et al. 


1993) of [2.2924,2.2977] pm. In AGNs this feature is diluted 


2009 Diamond-Stanic et al.||2009 LaMassa et 'al.|2Qiop. 

In summary, out of the full sample of 51 AGNs, we can 

determine the near-IR AGN fraction in 31 objects (for the 
others we give limits). Except for NGC 7496, all of these 
sources are also detected in the mid-IR and hard X-rays and 
all except four have been detected in [O IV]. This leads to 
a final sample size of 30 (27) for our multi-band AGN rela¬ 
tions (including [O IV]). Throughout this article we denote 
the AGN luminosity in the near-IR as derived from the di- 


in the nucleus due to the strong non-stellar continuum and 
we use this effect to measure the AGN contribution. 

The resulting EW maps are displayed in Eigs. [^and[^ 
for sources with and without detected dilution by non- 
stellar continuum respectively. Among the non-detections 
there are a few sources (NGC 613, NGC 3169, NGC 5643) 
which may show an AGN diluted core, but since the EW 
depression there is smaller than our adopted uncertainty 
we do not derive AGN luminosities from them. 

Given the EW map at each pixel x, we calculate the 
diluting radiation, i.e. pure continuum emission /dilute fhat 
does not show the CO absorption feature, using fdilute,x,y = 
EWint/EWa^^y where EWint is the intrinsic CO equivalent 
width (see below) of the galaxy outside the influence of the 
AGN light. The diluting radiation is then given as 


^ T he routine has alr eady been used for 
and Hicks et al. (2013). 


Davies et al, 


(2006) 
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Fig. 1: EW maps of sources with unambiguous nuclear dilution. We note that some of these sources (e.g. NGC 5135) 
show additional structure in the EW map. Article number, page 5 of 17 
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Fig. 2: EW maps of sources where nuclear dilution by AGN light cannot be safely established. 


-^dilute(^5 — -^total(^:^) -^stars(^:^) 

~ (1 1//dilute(^: 1/)) ' -^total(^5 P) (1) 

The fraction of non-stellar light in the near-IR is 

then calculated as the fraction of diluting radiation (over 
total) within an aperture A larger than the PSF: 


view) of these observations. The full-width at half maxi¬ 
mum (FWHM) of the point spread function (PSF) of these 
observations was 290 milli arcseconds (mas) and 181 mas, 
respectively. 

The adopted intrinsic EW per source and /ag^ are 
given in Tab. and its derivation is described in the next 
Section. 


^NIR _ 
JAGN — 


I A -^dilute (^: V) 
Ja -^totai 


( 2 ) 


We used an extraction aperture of 1" for all sources 
except for the Circinus galaxy and NGC 7469 where 0.5 
arcsec were used due to the smaller pixel scales (fields of 


The isotropic AGN luminosity in the near-IR, Ta^, 
is then simply given by /dilute * ^7rD‘^ where D is the dis¬ 
tance of the galaxy. These luminosities c an be compared 
to literature values, e.g. to Prieto et al. (2010) who have 
compiled a number of nuclear photometries. We find their 
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results for the K band photometry of the pure AGN source 
to be essentially identical to ours for most of the sources 


3.1.1. Determination of intrinsic CO equivalent width 

To determine the intrinsic, i.e. undiluted, equivalent width 
of the CO feature, we classify sources into three categories 
depending on the radial profile of their CO EW: 


irrelevant and the error is dominated by the photometric 
accuracy of our calibration which we also estimate to be 
about 10%. 

For sources where we cannot measure the intrinsic EW, 
we adopt a value of 11.1 A, the median of all undiluted 
sources, for the CO(2,0) equivalent width. This is slightly 
different from the value of 12.0 A that was suggested by 
( [20Q7| ). 


Davies et al 


1. Constant EW over the entire cube (i.e. undiluted by 
AGN light, EW = EWint), 

2. EW suppressed in the nucleus and more or less constant 
outside some radius (i.e. diluted by AGN light but in¬ 
trinsic EW is observed), 

3. constantly increasing EW (i.e. diluted by AGN light and 
field of view too small to observe EWint-)- 


The variety of EW profiles can be seen from the maps 
(Figs. [T||^and from the azimuthally-averaged radial pro¬ 
files (FigTW. One can see a clear bimodality in the nuclear 
value EWnuc- most sources are either highly diluted (EWnuc 
< 5 A) or essentially undiluted (EWnuc ^ EWint ~ 11.1). 
This is due to the fact that most of our sources have roughly 
the same stellar surface brightness (luminosity within our 
aperture), but we are probing a range of AGN luminosities 
(see Fig. [^. Depending on the AGN luminosity that we 
probe, the nucleus therefore appears either starlight domi¬ 
nated (undiluted CO feature) or AGN dominated (entirely 
diluted CO feature). The sources with intermediate values 
of EWnuc have AGN luminosities that are similar to the 
stellar luminosity. This transition occurs at ~ 10^^-^ 

erg/s which corresponds to Lboi ~ 10^^-^ erg/s. This hap¬ 
pens to be the luminosity below which the AGN torus is 
expec ted to disappear accord ing to some models ( [Elitznr 


2006 Honig & Beckert 2007). Unfortunately we can how- 


ever not constrain the torus disappearance with these mea¬ 
surements yet (we would need higher spatial resolution to 
probe lower AGN luminosities). 

For all sources except five, we can determine the in¬ 
trinsic value of EW from a fit to the radial profile. The 
median value is 10.7 A with a standard deviation of 1.6 
A (see Fig. [^. There is a slight difference between the di¬ 
luted sources (AGNs) and the undiluted sources with the 
former showing a median value oflO.3Tl.8A and the lat¬ 
ter 11.1 ± 1.1 A, but this difference is not significant given 
the relatively large errors. 

For each source where we can determine the intrinsic 
value, we use that value as a reference to calculate the frac¬ 
tion of AGN light. The accuracy of this method is limited 
by our ability to resolve small-scale EW variations which 
we estimate to be ^ 0.5 A. Thus we are effectively limited 
by spatial resolution, not by signal/noise. 

For stellar-light-dominated sources (low /ag^) this con¬ 
verts to an uncertainty in /ag^ of about 10%. For AGN- 
dominated sources, the uncertainty in the intrinsic EW is 


^ In NGG 1068, NGG 1566 and NGG 7582, however, our pho¬ 
tometries are significa ntly (about two mag nitudes or a factor 6) 
brighter than those of 


Prieto et al. 


(2010). This could indicate 


problems with the PSF correction in the very high spatial res¬ 


olution data from Prieto et al. (2010). Since our star-format ion 


corrected photometries are averaged over the central arcsecond, 
they should not be affected by PSF variations. 


8 



cn 

X 

2 


0 |_________ I ! |! |; _____ 

6 8 10 12 14 

intrinsic CO equivalent width [Angstrom] 


Fig. 3: Histogram of intrinsic equivalent width as deter¬ 
mined from Gaussian fits to the radial profile of those 
sources where a plateau level was seen. Black dashed line: 
median value of EW for all sources where the intrinsic EW 
can be measured (10.6 A); purple dashed line: median value 
for diluted sources only (10.3 A); orange dashed line: for 
non-diluted sources only (11.1 A). 


3.2. Spectral decomposition and near-IR color temperatures 

In addition to the starfit analysis, we also modeled the 
integrated nuclear spectrum (in an aperture of 1" for most 
galaxies, see above). We fitted a stellar template of a late- 
type star and a blackbody emitter (representing hot dust 
emission from the inner edge of the torus) to derive the color 
temperature of the continuum emission and an independent 
estimate of /ag^- An example of the spectral fit is shown in 
Fig.§ The resulting temperatures are mostly lower than 
the sublimation temperature of dust 1200 - 1800 K) 
and are given in Tab. Only a few galaxies show very 
blue nuclear spectra that are not consistent with thermal 
radiation from dust (no temperatures are given for these 
sources in the table). Most of these sources show no AGN 
light in our analysis of the CO equivalent widths. ESO 428- 
G14 and NGC 1097 are the only exceptions showing clear 
AGN signatures both in the near-IR, in X-rays and in the 
mid-IR, but a very blue continuum that is not compatible 
with thermal radiation from hot dust. 

The median temperature of Seyfert 1 (excluding li) is 
1292 ± 46 K, while Seyfert 2s show lower temperatures of 
887 ± 68 K. Sy li galaxies are in between at 1112 ± 47 K. 
The general trend of cooler K-band temperatures of Seyfert 
2s has been s een before, but mostly with somewhat lower 
temperatures. Alonso-Herrero et al. (1996) for example re- 
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Fig. 6. Fit to the nuclear (aperture 
1'') spectrum of NGC 1386. The ob¬ 
served spectrum is shown in blue, the con¬ 
tinuum (T=982 K) and stellar template 
(HD 176617) are shown as dashed and 
straight black lines respectively and the to¬ 
tal fit is shown in red. Prominent AGN 
emission lines and stellar absorption fea¬ 
tures are marked. In the inset we show 
our two-Gaussian fit to the Bry line of 
NGG 1386, demonstrating that the previ¬ 
ously found “broad” line is actually a su¬ 
perposition of two narrow lines. 



Fig. 4: EW averaged in annuli for all 51 sources. Only the 
good signal/noise part of the cube that was used to fit the 
EW profile is shown. Diluted sources are drawn in purple, 
undiluted in orange and sources where the intrinsic value 
could not be determined due to the limited field of view of 
the SINEONI data are shown in green. The black line de¬ 
notes NGC 5135 which has a non-axisymmetric EW profile. 
This graph shows sources observed in all three SINEONI 
pixel scales (field of view: 0.8, 3 and 8 arcsec). 


quire dust temperatures of 8QQ-12Q Q K to account for the 
AGNs in near-IR color-color plots; Riffel et al. (2009) re¬ 
port average temperatures of 1000 and 600 K for type 1 
and type 2 objects, respectively. The reason for their lower 
numbers is likely that they use lower spatial resolution data 
which are contaminated by cooler dust from regions off the 
nucleus. This is c onsistent with the lower AGN fractions 
(/nuc) reported by Alonso-Herrero et al. (1996). However, 
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Eig. 5: Distribution of AGN and stellar luminosities (surface 
brightnesses) within an aperture of 1 arc second; various 
levels of /agn indicated in the plot. 


the hotter dust temperatures we find are probably more re¬ 
alistic given the fact the the hottest dust is indeed at the 
sublimation radius as seen by near-infrare d interferometric 


and d ust reverberation observations (e.g. Kishimoto et al. 
2011). This innermost radius of dust is smaller than 1 m^ 
for the K-band brightest local Seyfert galaxies (Kishimoto 


et al. 2011), i.e. much smaller than our resolution limit of 
25 mas. However, all of these interferometrically studied 
sources show very high visibility in the K band, implying 
that the the continuum in the central few ten mas is dom¬ 
inated by hot dust emission, i.e. we are able to isolate this 
emission. 
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An additional result of this analysis is a completely inde¬ 
pendent estimate of the fraction of AGN light in the near-IR 
fAGN ~ same data and within the same aperture. We 

find very consistent results for most sources (Fig. 0. Only 
in one source (NGC 1386), the results deviate substantially 
if AGN = 12 % in the radial decomposition, 25 % in the 
spectral decomposition). 

It is perhaps unexpected that these two methods work 
about equally well when one may have expected the ra¬ 
dial decomposition to be a more solid analysis. The latter 
relies only on the spatial invariance of the CO EW over 
the central 100 - 200 pc (which, empirically, turns out to 
be the case) and not on the value of CO EW itself. It is 
helped by looking at the 2.3 /rm CO bandhead, the depth 
of which does not vary much with star formation history for 
an evolving s tellar population w ith only a few specific ex¬ 
ceptions (e.g. Davies et al.|2007 ). The radial decomposition 
only uses the CO EW ana is thus insensitive to variations 
in the spectral slope, i.e. extinction. 

The spectral decomposition on the other hand has the 
observational advantage of not requiring spatial informa¬ 
tion at all, but relies on the invariance of both the spectral 
slope and the actual value of the CO EW (i.e. the stellar 
template chosen). It is sensitive to extinction, which can 
make the spectral slope redder (although that is mitigated 
by the fact we are looking in the K-band). This method 
will over-estimate the AGN fraction when there is signifi¬ 
cant reddening to the stellar population. This may be the 
reason why the /ag^ values from the spectral decomposi¬ 
tion tend to lie slightly above th e va lue from the radial de¬ 
composition in Eig.[^ In Section [4^ we discuss the relation 
between reddening of the AGN light and color temperature 
of the AGN heated dust in the context of a simple AGN 
obscuration model. 


In conclusion it seems that both methods give similarly 
good results for the decomposition of stellar and AGN light. 
Given the fact that the dust temperatures derived from 
spectral decompositions are often higher than the dust sub¬ 
limation temperature, however, we would still recommend 
to be cautious when only a spectral decomposition is used. 
If possible, it is advisable to use both the spectral and radial 
decomposition (on lEU data) to estimate /ag^- 


3.3. AGN classifications 


AGN fractions in the near-infrared (/ag^) higher for 
optical type Is than for optical type 2s, as expected, but this 
effect is entirely due to a luminosity mismatch between the 
two sub-samples as shown in Eig. There seems to be no 
difference between /ag^ of Wpo I and that of type 2 sources 
at any specific luminosity. The fact that /ag^ can differ 
substantially for AGNs in our luminosity range means that 
one has to be very careful when constructing AGN SEDs 
from large-aperture photometry data in the near-infrared. 

Here we classify AGNs as type 2 if they belong to the 
NED classes ”Sy2" or "Sylh” (Seyferts with hidden broad 
line regions that can be seen in polarized light). Classes 
Syl, Syl.5, Syl.8, Syl.9, Syln and Syli will be put into 
the type I bin - but we will mark Syli specially to show 
how important it is to correctly classify these intermediate 
sources as type I in infrared studies. 


Eor NGC 1386 the detection of a somewhat broad 


800 km/s) Brq line has been reported by Reunanen et al. 



Eig. 7: Comparison of /ag^ values (within I arcsecond for 
most sources) as derived from the radial profile of the CO 
EW and from the spectral decomposition. 


(2002) and the source is therefore listed as an infrared type 
1 Seyfert galaxy in NED (although most Syli have line 
widths ^ 1000 km/s). With our higher spatial resolution 
and better S/N data, we cannot confirm this broad line, but 
instead find that the Bry line is composed of two narrow 
(^ 300 km/s) components that are offset by 355 km/s (see 
Fig. 0. We will therefore treat this galaxy as a Seyfert 2 
galaxy in our analysis. 
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Eig. 8: Near-IR AGN fraction /ag^ ^.s a function of near-IR 
AGN luminosity for the various classes of AGNs. The ap¬ 
parent correlation between /ag^ and Tagn (NIR) is caused 
by the fact that the stellar luminosity is essentially constant 
within our aperture in all objects. 
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Fig. 9: AGN luminosity relations. From top to bottom, left to right: (1) mid-IR - near-IR, (2) hard X-rays - near-IR, 
(3) mid-IR - hard X-rays, (4) mid-IR - [OIV], (5) hard X-rays - [OIV]. Luminosities are given in log (L/[erg/s]). The 
straight lines indicate the best fit, the dashed lines the ~ la uncertainties and the dotted lines the best fit with slope 
fixed at 1; see text for details of the fitting procedure. 


4. Discussion 

4.1. AGN luminosity relations 


Linear fits and correlation analysis Having compiled a large 
sample of nuclear AGN luminosities in a number of bands, 
we will now discuss their relations (Fig.[^. For each relation 
(in log-log space), we generate 10^ re-sampled data points 
according to their two-dimensional (assumed Gaussian) er¬ 
ror distribution and compute linear fits using a simple 
metric. We use the median of the resulting distributions to 
represent the best fit relation and cast our lower and upper 
uncertainties in terms of the 16^^ and 84^^ percentiles of 
the distributions. The results are given in Tab. 


NGC 1068 has been excluded from all fits involving X- 


ray m easurements since it is heavily absorbed (Bauer et al. 


2014). This affects even the very high energy band used by 


Swift/BAT as can easily be seen from its large offset from 
the respective relations. 


4.1.1. Near-IR luminosity relations 

We find that a strong correlation exists between the nu¬ 
clear near-IR and mid-infrared luminosities of local AGNs 
(Fig. top left) as well as between the near-IR and hard 
X-ray luminosities (Fig.l^ top right). The best fit is slightly 
different from a linear rSation which is indicated as a dot¬ 
ted line in the fits. However, we consider the differences 
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marginal and note that there is no indication in the litera¬ 
ture of a non-linear behavior in the luminosity r ange that 


W 

th< 


we cover (Lutz et ah 2004 Gandhi et ah 2009). Only at 


higher luminosities, the X-ray obscured fraction decreases 
(Merloni et ah 2014) which can be very well explaine d in the 


context of t he receaing (in fact: sublimating) torus (Davies 
|et al.|[201^ . This may change the slope (especially for the 
mid-IR - X-ray relation, see below) at X-ray luminosities 
> 10^^ erg/s. 

Therefore, and since the scatter is dominated by intrin¬ 
sic variations rather than observational errors, we compute 
the two-sided Kolmogorov-Smirnov statistic of the ratio of 
the luminosities to determine the significance of the offset 
between type 1 and type 2 AGNs. Both near-IR relations 
show a significant offset between type 1 and type 2 sources. 

^AGN (^agn) = 1 sources are about 7 (10) 

times brighter in the near-IR than type 2 sources with very 
high significance (see Tab.|^. 

I t is important to note that infrared Seyfert 1 galaxies 
(e.g. |Lutz et al.||2002|), where broad lines are visible only 
in infrared bands, lie in the same locus as optical type 1 
sources and should therefore be treated as type 1 objects 
in infrared studies. This is perhaps not surprising since our 
analysis is based upon infrared emission as well, but it is 
an often overlooked issue. The K-S distance and probability 
indicate a better separation and higher significance of the 
separation if type li objects are treated as type 1 AGNs 
(see Tab. [^. We have therefore marked these sources as 
type 1 objects (blue) with a red box around the symbol. In 
our type 1 sample, five out of 16 sources belong to the li 
class. 

Ou r resu lts can be compared to the finding from [Honigl 


et al. (2011) who estimated the anisotropy in the infrared 


emission in a sample of powerful quasars and radio galaxies 
and find an anisotropy factor of about eight at 2.3 jam which 
is consistent with the offset we find in the near-IR - mid-IR 
relation for our sample of much lower luminosity AGNs. 

4.1.2. Isotropic AGN luminosity relations 

From the same sample, we also construct the relations be¬ 
tween mid-IR, hard X-ray and [O IV] luminosities (Fig. 
bottom row), all of which are supposed t o be indicators 
of th e isotropic luminosity of AGNs (e.g. LaMassa et al.| 
2010). In our compilation we indeed do not see significant 
differences between type 1 and type 2 sources in these wave¬ 
bands. The best fits to the relations indicate that some 
of them are highly non-linear (see Tab. [^. Due to the 
larger scatter, however, we do not further discuss these non- 
linearities. In the plots we show the best fit for an assumed 
linear relation. 

For the well-known mid-IR - hard X-rays relation there 
is a marginal indication (on a 10% significance level) that 
type 2 sources are brighter by a factor of about 2 in 
at 42.5 than type 1 sources. This indicates that the 

mid-IR emission from the AGN “torus” is indeed slightly 
anisotro pic as suggested by bo th smooth and clumpy torus 
models. |Nenkova et al. (2008), for example, find a factor 2 
incr ease in flux between edge-on and face-on clumpy tori 
and Snyder et al. (2013) find the least obscured sources to 
be more than a factor of ten more luminous at 12 /im than 
the most obscured sources. In most studie s of this relation, 
however, no difference has been seen (e.g. Lutz et aT]|2004 


Gandhi et al.|2009 Asmus et al.|2011 ). We believe that the 
marginal difference we see arises due to two improvements 
over previous studies: 

1. Absorption corrected 2-10 keV luminosities are a less re¬ 
liable indicator of the hard X-ray radiation of an AGN 
than the Swift/BAT fluxes from the 70-month data re¬ 
lease that we mostly use (for 24 out of 30 sources). The 
Swift/BAT fluxes are essentially unaffected by absorp¬ 
tion (the attenuation is less than 2 for log{N h/ cm^) < 
24.3). Due to the long-term nature of the BAT survey, 
variability effects are also reduced. 

2. The small anisotropy in the hard X-ray - mid-IR rela¬ 
tion is washed out if type li sources are not classified as 
type 1 (“face-on”) sources. 


Most recently, [Mateos et al. (2015) studied the rela¬ 
tionship between the 6 /im and 2-10 keV luminosities of 
a sample of relatively powerful AGNs and also found an 
anisotropy in the sense that, at a given 2-10 keV luminos¬ 
ity, type 2 AGNs are ^ 1.3-2 times fainter at 6 /im. This 
indicates that the 6 jam emitting region is more closely re¬ 
lated to the 12 /im emitting region (which has a similar 
anisotropy) than to the hottest dust, emitting at 2.3 /im, 
which has a much larger anisotropy of a factor of about 
7. This is e ncouraging for the upcoming VLTI instrument 
MATISSE ( Lopez et al.pOM ) which will open the 3-5 /im 
window for imaging interferometry. The small anisotropy 
at 6 /im suggests that at the wavebands of MATISSE, one 
will be able to resolve the relatively large warm body of the 
torus and not the innermost hot dust (which is essentially 
unresolved with the currently available baselines). 

There is a slight offset between type 1 and type 2 sources 
in the [O IV] - mid-IR relation (Eig.[^ bottom center) that 
would also indic ate that the mid-I R emission is anisotropic 
(as also seen by [Yang et al.j (|2015|)). However, this offset is 
not significant. There is nearly no difference between type 
1 and type 2 objects in the [O IV] - hard X-ray relation 
(Fig.§ bottom right). 

4.2. Two component model 

If both type I and type 2 sources are drawn from the same 
population - the central premise of unification models - 
then the clear offset of the type I and type 2 sources in the 
near-IR - mid-IR luminosity relation has to be explained 
simply by inclination, i.e. different amounts of obscuration 
of the central source. 

In order to test the effect of obscuration, we plot the 
near-IR color temperatures that we derive from our spec¬ 
tral fitting against the offset from the near-IR - mid-IR 
luminosity relation (Eig. 10). Here we have converted this 
offset to an astronomicalmolor” so that it can be easier 
compared to other observations. We refer to the near-to- 
mid-infrared dust color as “K-N” although the magnitudes 
are actually not drawn from band-averaged fluxes, but from 
2.3/im and I2.0/im monochromatic fluxes, respectively. 

The first thing to note in our temperature-color plot is 
that AGNs of all optical types have very red near-to-mid- 
infrared colors: K — N > 4. As expected, type Is are less 
red and hotter than type 2s, but both have near-infrared 
temperatures that are indicative of hot dust. Again we see 
that infrared type I Seyfert galaxies are closer to optical 
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K-N 


Fig. 10. Near-IR color temperature vs. 
NIR-MIR color. Red infrared colors are 
to the right. The arrows are extinction 
vectors from our two blackbody model 
(crosses are at Ay^^ = 0,5,10,... mag). 
Solid arrow: Thot = 1500 K, Twarm = 300 
K, /area = 1000. The dotted vectors are for 
a third (left) and three times (right) /area- 
The curled (dash-dotted) vector is with a 
scattering fraction of 10% (for details see 
text). 


type 1 than type 2 galaxies in this plot with temperatures 
intermediate between the two classes. 

To explain both the hot near-IR temperatures and the 
red near-to-mid-infrared color, we require a hot and a 
warm component and only the hot one can be signifi¬ 
cantly inclination dependent (obscured) because of the near 
isotropy of the mid-IR radiation. Such a two-component 
model is also observationally motivated as AGN tori are 
known to consist of tw o spatially distinct com ponents of 
hot and warm dust (e.g. Bnrtscher et al.|[2QQ9 ). 

The intensity in this model is then given by 


h = 


rhot 


+ /a 


rwarm 


(3) 


where /area is the area ratio between the warm and the hot 
emitter and the two components are given by 


rhot _ 

\ — 


Bx{X, Thot) 


exp(— rv(A) • 


V /1-09), 

Twarm)- 


(4) 

(5) 

( 6 ) 


Here Bx is the Planck function for the spectral radiance 
in wavelength space, Thot and Twarm are the temperatures 
of the respective components and ry(A) = ry/rv is the 
optical depth of the dust screen at a given wavelength. We 
interpolate it loga r ithmi cally from the extinction curve of 
Schartmann et al. ( |20Q5[ ) and give it relative to the optical 
depth at V band (5500 A) where it is normalized to 1. 
is the screen extinction towards the hot component at V 
band in magnitudes. We do not require an extinction to 
the warm component. 

The temperatures of our simple two-component “torus” 
model are set to Thot = 1500 K and Twarm = 300 K. The 
hot temperature is chosen according to the hottest tem¬ 
peratures that we measure from the near-IR spectra. The 
warm temperature is not very well defined by this analysis, 
but since the maximum of the re-emitted AGN radiation is 
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at ~ 10 /im, the temperature of the warm co mponent must 
be at about 300 K (see e.g. Figs. 3 and 4 of Prieto et al. 

( Mq| )). 

The K-N color of the type 1 galaxies in our temperature- 
color plot is set by /area ~ 1000. This value can be com¬ 
pared to interferometrically measured radii for the warm 
and hot dust, the ratio of which is 10 ... 30 (see Fig. 36 of 


Bnrtscher et al. ( 2013[ )), i.e. /area = 100... 1000, in reason- 
able agreement with our temperature-color plot. Fig. 
where the horizontal range covered by the dotted lines cor¬ 
responds to /area = 300 ... 3000. The modeled range of /area 
is about three times larger than what one may expect from 
the interferometric observations. This difference cannot be 
explained by extinction to the mid-IR component, since an 
obscuration of about 2 mag (at 12 jam) would have the 
same effect as reducing f^rea by a factor of 3. However, a 
slight variation of Twarm by a factor of 3^/^ « 1.3 suffices 
to explain the small discrepancy. 

The only other free parameter of our model is the ex¬ 
tinction to the hot component which may be related 
to the inclination of the “torus” axis out of our line of sight. 
This parameter takes values of 0... 5 mag for optical type 
1 Seyferts, 5... 15 mag for type li Seyferts and 15... 35 
mag for type 2 S eyferts. Since ry 13 • r 2 . 3 ^m (Schart¬ 


mann et al. 2005), the extinction for type li Seyferts can 
also be expressed as ^ 2 ° 3 Vm = 0.4 ... 1.2 mag. This is con¬ 
sistent with their classification which is based on the fact 
the broad lines are visible in the infrared. Probably these 
sources host an intrinsic type 1 nucleus but are slightly ob¬ 
scured by foreground dust in their host galaxy. 


Even this simple model is degenerate: a variation in 
Twarm bas nearly the same effect as a variation in /area and 
instead of obscuring the hot component by small amounts 
of dust (i.e. moderate inclinations?), we could also attribute 
the range of Seyfert 1 temperatures to intrinsic variations. 
But the point of this model is that with only two free pa¬ 
rameters (/area: ^v°^) as wcll as two obscrvatioually de¬ 
termined parameters (Thot^Twarm) we can explain almost 
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the entire range of near-IR temperatures and near-to-mid- 
infrared colors in our sample of local Seyfert galaxies. 

There are four outliers that cannot be explained by a 
variable amount of extinction of the hot component and 
small variations of /area- On the blue side there is one 
source, ESO 548-G81 which may just have a lower frac¬ 
tion of warm to hot dust compared to the bulk of Seyfert 1 
sources. On the red side there are three outliers: NGC 1386, 
IC 5063 and the Circinus galaxy. In the reddest of all, the 
Circinus galaxy with K-N 10, the hot component seems to 
be intrinsically much weaker than in the other sources. This 


has been seen in high-resolution SED models (Prieto et al. 


2004) and recently been confirmed by a det ailed analysis of 


the highest resolution interferometric data (Tristram et al. 


2014). Although the Circinus galaxy is so nearby that the 


sublimation radius is nearly resolved, there is no evidence 
for hot dust in the interferometric data. This is in contrast 
to NGC 1068, where a similar resolution is reached in in¬ 
frared interferom etric observation and a hot component is 
clearly detected (Lopez-Gonzaga et al. 2014) - consistent 
with the “bluer” color of NGC 1068 compared to the Circi¬ 
nus galaxy. 

The other two red outliers, NGC 1386 and IC 5063, are 
harder to explain as they are both red and hoi|^ There is 
no extinction vector that would connect the locus of these 
sources with a plausible “starting” position among the type 
1 sources, but there are a few possible ways of increasing 
the color temperature in the K band while keeping the red 
colors. 

A significant contribution from a non-thermal source 
with a flat spectrum oc would increase the modeled 
temperature, but most Seyfert galaxies are not dominated 
by non-thermal radiation in the infrared, as evidenced for 
example by the fact that the infrar ed emission is resolved 
in most objects with interferometry ( Burtscher et al.|2Q13 ). 

Alternatively, if some part of the hot dust radiation 
could escape unreddened, it would also affect the K band 
color in the desired way, i.e. if we replace by 


+ fsct ■ Ix{X,Thot), 


(7) 


where /scat is the fraction of hot blackbody emission 
that is scattered into our line of sight without obscuration. 

One possibility to achieve this would be via elastic scat¬ 
tering off free electrons in the narrow-line region. This is 
perhaps attractive since IC 5063 is observed to have a 


hidden broad-line region as seen in polarized light (Inglis 


et al. 1993). This indicates that a scatterer exists in this 


source. The scattered radiation should be measurable via 
a high degree of polarization also in the near-infrared as 
the Thomson cross-section for scattering off free electrons 
is wavelength-independent. Eor an extinction vector includ¬ 
ing scattering to match the red and hot outliers in Eig. 
(the curled dot-dashed line), we require a scattering frac¬ 
tion of 10 % which corresponds to an intrinsic polarization 
of about 60 % in IC 5063 since the direct radiation would 


^ NGC 1386 should perhaps be treated with caution as it is the 
only source with a factor two difference in /agn between our two 
analyses. We have plotted the K-N color derived from the radial 
profile decomposition, but note that the source would move to 
the left by about 0.8 mag if we had taken the value derived from 
the spectral decomposition. We have, however, no reason not to 
trust the radial decomposition in this source. 


be attenuated by A 2 . 3 ^m ~ 3. Eor this so urce, very detailed 
polarization observations are available by [Lopez-Rodriguez| 


et al. (2013) who measure a degree of polarization of 8 % at 
2.3 jam. Using our own /agn aperture correc¬ 

tion factor of 1.6, this corresponds to an intrinsic polariza¬ 
tion of 23 %. Scattering therefore seems t o be an unlikely 


explanation. Lopez-Rodriguez et al. (2013) also considered 
synchrotron radiation as the origin of the polarized light 
in IC 5063, but since the source shows no variability, it is 
unlikely to host a strong non-thermal source. 

Another way of “by-passing” the obscuring torus would 
be to give up the assumption of uniform screen absorption 
of the hot component. This could be through “holes” in the 
clumpy torus or because only a part of the region hosting 
the hottest dust is covered by colder torus clumps or be¬ 
cause we are seeing the source at a grazing inclination angle. 
This may imply that these sources are good candidates for 
“changing-look” AGNs which are in the proces s of changing 


betwe en obscured and unobscured classes (e.g. Risaliti et al. 


2010). However, the red outliers in the temperature-color 


plot are all (moderately) Compton-thick sources, where one 
may not expect holes in the torus. In conclusion, most 
the sources are explained well with this very simple two- 
component model and the few outliers probably represent 
the intrinsic diversity in AGN nuclear dust structures. 


5. Conclusions 

Eor an archival sample of 51 local Seyfert galaxies, we ana¬ 
lyze high-resolution near-IR integral field spectroscopic ob¬ 
servations, mostly from VLT/SINEONI observations. We 
use these data to produce equivalent width maps of the 
stellar CO A2.29/im absorption feature. Erom the dilu¬ 
tion of this feature, we infer the AGN fraction and the 
AGN luminosity in the near-IR i.e. hot dust lu¬ 

minosity. We find a significant fraction of nonstellar light 
in 31 sources. Erom the literature, we additionally collect 
measurements of the nuclear mid-IR radiation from high- 
resolution ground-based observations, hard X-ray fluxes 
from the Swift/BAT all-sky survey and other observations 
as well as [OIV]A25.89/im fluxes from Spitzer observations. 
These measurements are considered isotropic indicators for 
the AGN luminosity. Equipped with this collection of high 
resolution data, our findings are: 

— The near-IR AGN fraction /agn (within 1 arcsec) varies 
quite significantly among AGNs. In our sample, type 1 
Seyfert galaxies have higher /agn Itian type 2s only 
because of a luminosity mismatch between the two sub¬ 
samples. AGNs of both types can take on fractions from 
near 0 to 100 %. It is therefore important to take into 
account contamination by stellar light even when ob¬ 
serving with small apertures. 

— There is an apparent bimodality in the nuclear equiva¬ 
lent width of the stellar CO feature in nearby galaxies. 
Inactive galaxies and weak AGNs show an almost undi¬ 
luted feature while powerful AGNs are strongly diluting 
the stellar light. The transition occurs at the luminosity 
where the torus is expected to vanish. However this is 
a coincidence that is caused by the fact that the stellar 
surface brightness is nearly constant for our sample of 
galaxies over a wide range of AGN luminosities. 
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— Using our estimate of the near-IR AGN luminosity with 

high-resolution mid-IR and hard X-ray measurements, 
we find a strong correlation with significant offsets be¬ 
tween optical (and infrared) type 1 AGNs and type 2 
AGNs. They are typically 7 (10) times brighter in the 
NIR at a typical luminosity of log =42.5 (log 

^AGN = 42.5). 

— We also find a marginally significant offset between type 
1 and type 2 sources in the hard X-rays-mid-IR re¬ 
lation in the sense that type 1 sources are about two 
times brighter at a typical hard X-ray luminosity and 
attribute this to anisotropy in the mid-IR emission of 
the torus. The relation between [O IV] luminosity and 
hard X-ray luminosity, on the other hand, shows almost 
no offset between the two classes of AGNs. This suggests 
that [O IV] and hard X-rays are equally isotropic indi¬ 
cators of the AGN luminosity and the anisotropy in the 
other relations arises in the near- and mid-IR emitting 
regions. 

— We invoke a simple two component model that consists 
of a warm (T^arm ~ 300 K) and a hot (Thot ~ 1500 K) 
blackbody and an area scaling factor between them. We 
find that the scaling factor (^ 1000) is consistent with 
the expected area ratio given the interferometrically ob¬ 
served ratio of mid-IR to near-IR radii. Almost all our 
sources can then be explained by simply adding an ab¬ 
sorbing screen to the hot component. An Ay of 5-15 
mag is required for type li sources (corresponding to ^ 

1 mag in the near-IR) and 15-35 mag is found for type 

2 sources; pure type 1 sources cover an Ay range of 0-5 
mag. There are only four outliers to this simple scenario. 
Three of them are Seyfert 2 galaxies with a very red K-N 
color but hot near-IR color temperature. This can possi¬ 
bly be explained by scattering, non-thermal emission or 
variations from a screen absorber, but detailed studies 
are needed to clarify this. 

— We find that infrared type 1 Seyfert galaxies, where 
broad lines are only visible in infrared wavebands (such 
as Bry), are more similar to type 1 galaxies also in other 
infrared indicators (e.g. /ag^ and K band color temper¬ 
ature) and suggest to treat these Seyfert li galaxies as 
type 1 sources. A mis-classification of these sources may 
otherwise complicate results from infrared studies. 
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Table 2: EWint- intrinsic equivalent width of the CO feature. For sources where it cannot be derived from the data due to 
the limited field of view, the median value for non-diluted cubes (11.1 A) is used. Luminosities are given in log(L/[erg/s]). 
Mid-IR, X-ray and [O IV] luminosities are only listed for sources where we can measure the NIR luminosity of the AGN. 
The K magnitude is extracted within 1 arcsecond, e xcept for the two sou rces marked with ^ where t he a perture is 0.5 

( 2Q08D, b: Burtscher et al. J2013[ ), c: [Asmus et ah 
^rtner et al.||2Q13|), absorption-corrected 


arcseconds. References for mid-IR measurements: a: Raban et al. 


(2014). Refe rences for X-ray lumin osi ties: BAT: 70-mQnth B AT survey 
Levenson et al.|(2Q06), e: Gandhi et al 


3aum; 


2TTQlceV: d: 
Diamond-Stanic et al 


(|2QQ9), h: [Weaver et al. "pOToj, i: [Tommasin et al.|(|2QT^ 


(20091), hlAsmus et al.l(|201l|). References for [O IV] luminosities: g: 
r|Tommasinetal.|(|2010|) 


Source id 

K mag 

[A] 


T[K] 

log (^a™) 

log 

log (^agn) 

log(iS) 

Gircinus ^ 

11.11 

(11.1) 

77 

691 ±2 

40.73 

42.73^ 

41. 

40.16^ 

ES0428 

11.75 

10.5 

14 

- 

41.27 

42.31^ 

43.20'^ 

- 

ES0548 

10.30 

7.2 

87 

1236 ± 3 

43.39 

42.95^ 

43.26®^'^ 

- 

IC1459 

11.54 

9.6 

< 10 

- 

<41.33 

- 

- 

- 

IC5063 

12.89 

7.8 

56 

947 ± 11 

41.87 

43.63^ 

43.18^^'^ 

41.39^ 

M87 

12.56 

12.8 

66 

1423 ± 8 

41.30 

41.27^ 

41.43/ 

- 

MCG523 

10.43 

8.0 

86 

1207 ±5 

43.00 

43.53^ 

43.58®^"^ 

40.72'* 

NGC289 

13.17 

11.8 

< 10 

1006 ± 61 

<40.44 

- 

- 

- 

NGC613 

12.94 

10.7 

< 10 

- 

<40.46 

41.53^ 

- 

- 

NGC676 

12.16 

11.2 

< 10 

- 

<40.75 

- 

- 

- 

NGG1052 

11.43 

11.2 

24 

1239 ± 14 

41.45 

42.22^ 

42.15®^'^ 

39.06^ 

NGG1068 

7.90 

7.5 

91 

723 ± 1 

43.16 

44.00^ 

41 94BAT 

41.67^ 

NGG1097 

12.09 

10.7 

25 

- 

41.08 

41.17^ 

41.50® 

39.3T 

NGG1365 

9.73 

6.4 

94 

1282 ± 4 

42.63 

42.56^ 

42.39®^'^ 

40.78^ 

NGG1386 

11.85 

10.4 

12 

982 ± 12 

40.91 

42.48^ 

42.80/ 

40.54^ 

NGC1566 

11.38 

9.1 

56 

1479 ± 8 

41.42 

41.42^ 

41.58^^'^ 

39.20^ 

NGC2110 

10.89 

9.5 

80 

967 ±2 

42.69 

43.08^ 

43.69^^'^ 

40.84^ 

NGG2911 

12.79 

11.1 

< 10 

- 

<41.35 

- 

- 

- 

NGG2974 

11.80 

11.0 

< 10 

1106 ±32 

<41.10 

- 

- 

- 

NGC2992 

11.57 

12.6 

61 

1111 ±7 

42.20 

42.76^ 

42.51^^'^ 

41.IF 

NGG3081 

12.83 

12.8 

36 

977 ± 10 

41.32 

42.53^ 

42.85®^"^ 

40.92^ 

NGG3169 

11.46 

10.2 

< 10 

1600 ± 40 

<41.00 

40.93^ 

- 

- 

NGC3227 

10.53 

9.7 

67 

1543 ± 15 

42.29 

42.62^ 

42_76BAt 

40.47^ 

NGC3281 

12.06 

9.2 

82 

732 ±3 

42.51 

43.36^ 

43.39^^'^ 

41.60^ 

NGC3312 

12.44 

11.4 

< 10 

518 ±15 

<41.36 

41.68^ 

- 

- 

NGG3627 

11.56 

11.3 

< 10 

1433 ± 39 

<40.06 

40.23^ 

- 

- 

NGG3783 

10.00 

(11.1) 

99 

1244 ± 3 

43.39 

43.67^ 

43.70®^'^ 

40.88^ 

NGC4051 

11.14 

8.8 

64 

1305 ± 4 

41.64 

42.39^ 

41 92BAT 

39.75^ 

NGG4261 

12.99 

9.6 

< 10 

1455 ± 22 

<40.78 

41.53^ 

- 

- 

NGG4303 

12.39 

13.8 

< 10 

702 ± 22 

<40.44 

40.62^ 

- 

- 

NGC4388 

12.09 

10.7 

77 

887 ±4 

41.75 

42.41^ 

43.18^^^ 

41.15^ 

NGG4438 

11.63 

11.0 

< 10 

1185 ±33 

<40.68 

40.77^ 

- 

- 

NGG4472 

12.88 

9.0 

< 10 

1603 ± 28 

<40.35 

- 

- 

- 

NGC4501 

11.75 

10.9 

< 10 

- 

<40.79 

40.53^ 

- 

- 

NGG4569 

10.92 

11.2 

< 10 

- 

<40.86 

- 

- 

- 

NGG4579 

11.34 

10.3 

33 

1374 ± 12 

41.67 

42.00^ 

41.80^ 

39.18^ 

NGG4593 

11.02 

(11.1) 

87 

1292 ± 4 

42.42 

42.72^ 

42 87BAT 

40.05^ 

NGC4762 

12.09 

11.4 

< 10 

- 

<40.65 

- 

- 

- 

NGC5128 

10.20 

10.3 

81 

796 ±1 

41.03 

41.84^ 

42.39^^"^ 

39.24^ 

NGG5135 

13.03 

13.0 

54 

1033 ± 9 

42.12 

43.15^ 

43.70® 

41.40^ 

NGG5506 

9.25 

(11.1) 

99 

1223 ± 5 

43.09 

43.22^ 

43.22^^"^ 

41.18^ 

NGC5643 

12.08 

10.1 

< 10 

900 ± 30 

<40.66 

42.34^ 

41 79BAT 

- 

NGG6300 

12.05 

9.9 

53 

677 ±2 

41.21 

42.48^ 

42 34BAT 

39.82^ 

NGG6814 

11.98 

11.7 

58 

1049 ± 4 

41.73 

42.18^ 

42.68®^'^ 

40.12^ 

NGC7130 

13.06 

12.7 

< 10 

1061 ± 56 

<41.42 

43.11^ 

42.89^^"^ 

- 

NGG7135 

12.20 

10.5 

< 10 

- 

<41.24 

- 

- 

- 

NGC7172 

11.22 

9.6 

70 

822 ±3 

42.46 

42.80^ 

43.37^^'^ 

40.83^ 

NGC7469 * 

10.80 

(11.1) 

88 

1355 ± 6 

43.17 

43.78^ 

43.411^^'^ 

41.15^ 

NGC7496 

12.62 

13.1 

0 

667 ± 53 

38.75 

- 

- 

- 

NGG7582 

9.63 

10.5 

89 

1082 ± 2 

42.78 

42.95"^ 

42.63^^^ 

41.07^ 

NGC7743 

12.15 

12.1 

< 10 

1563 ± 82 

<40.75 

- 

- 
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Table 3: Fit results for various AGN luminosity relations. The intercept and slope for e.g. the NIR-MIR relation is defined 
as: logl/^^ = intercept + slope x logL^Q^. For each relation, the K-S test is computed from the ratio of luminosities 
for type 1 vs. type 2 AGNs. 


relation 

intercept type 1 

slope type 1 

intercept type 2 

slope type 2 

K-S distance 

K-S prob 



type li treated as type 1 




NIR-MIR 

5.287t{iS39 

0.870t°;°29 

1.844tP24I 

0 Q91+0-029 
o.aoi_Q q29 

0.854 

2.5e-05 

NIR-X-rays 

4-748t{;^« 

n CQ9+h046 

U.obz_Q Q4Q 


0.806t°;jf9 

0.875 

2.5e-05 

X-rays-MIR 

5.308t|367 

n C77+h054 

9111 q+4-468 

Zl.l±0_4 4Q4 


0.449 

l.Oe-01 

OIV-MIR 

o iro7-|-l-299 
o.00/_]^ 289 

0 747+0-030 
-0.030 

9 i^Q7+i-030 

Z.oy/_! 655 

1 ni 9+0-038 

i.UiZ_Q Q3g 

0.292 

5.9e-01 

OlV-X-ray 

10.673ti;14^ 

O.Oy /o.o27 

0.968t®;^“ 

O.yZD 0.i36 

0.286 

6.5e-01 



type li treated as type 2 




NIR-MIR 

4.624tPJ8 

0.886^^29 

-5.062tJi44 

1.097tni? 

0.610 

7.3e-03 

NIR-X-rays 

0.586t2;438 

0.980^^57 

5.467^^^39 

0.848t°;J^2 

0.722 

9.2e-04 

X-rays-MIR 


0.8481^^7 

17.61V3;471 

0.592tnri 

0.375 

2.5e-01 

OIV-MIR 

10.115tJio2 

0.708tn!J 

2.813tt722 

0.886^^40 

0.340 

4.3e-01 

OlV-X-ray 

6.932tPJ? 

n 7Q9+h036 

U. 1 OZ_Q Q3g 

13.078t|3i^ 

0.645^^6? 

0.422 

2.0e-01 
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